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Abstract 

We have investigated effects of an external magnetic field in the impurity 
Anderson model with a finite on-site Coulomb repulsion U. Large Nf expan- 
sion is employed in the slave boson representation, by taking into account / , 
/ , and / 2 subspaces. To evaluate the vertex function for the "empty state 
boson" self-energy, we have devised two approximations which reduce much 
computational efforts without losing general features of the model. It is found 
that the Kondo temperature is reduced by the presence of a magnetic field, 
and that at low field and at low temperature, the field dependence of both 
the Kondo temperature and the impurity magnetization exhibits a scaling 
behavior with high accuracy. Further, some interesting features are found 
in the field dependence of the impurity magnetization at finite temperature, 
the physical implications of which are discussed in terms of the renormalized 
Kondo temperature. 

PACS numbers: 71.10,+x, 71.27.+a, 71.28.+d, 75.20.Hr 
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I. INTRODUCTION 



Heavy fermion and valence fluctuation phenomena (see e.g., observed in Ce and Yb 
compounds have been often described by the impurity Anderson model ||. Due to a large 
on-site Coulomb repulsion U between /-electrons, these compounds are considered as typical 
examples of strongly correlated systems to which the conventional canonical perturbation 
technique cannot be applied. A systematic large- Nf (Nf the spin and orbital degeneracy) 
expansion method has been used for the assumed infinite-?/ model to describe dynamical 
as well as thermodynamic properties The infinite-?/ treatment allows one to consider 

only f° and f 1 subspaces of /-electrons, and thus makes the problem quite simple. However, 
this scheme cannot be justified for realistic systems because U is not an infinity (5~ 6eV 
typically) . Especially in Ce compounds, f 2 configurations are even energetically comparable 
to the /° configuration. 

In the "spin fluctuation" limit (or Kondo limit), the finite-?/ effect is manifested in vari- 
ous physical properties through the renormalized Kondo temperature Ta- The renormaliza- 
tion of Ta arises from the change of the exchange coupling constant J between conduction 
electron and /-electron spins due to a finite U, in a sense of the Schrieffer- Wolff transfor- 
mation . While this effect was recognized in the exact Bethe ansatz approach and also 
in the numerical renormalization group approach || for Nf = 2, more appropriate descrip- 
tions of the systems with Nf > 2 were provided by the large- Nf treatments including the 
f 2 - subspace in the variational approach fl||| and in the perturbational approach flOHT5 



Using the generalized slave boson technique, Schiller and Zevin |TJ| have shown that the 
lowest order perturbation expansion reproduces the exact Nf — ► oo variational result at zero 
temperature (T = 0) |§. 

In this paper, we have investigated effects of an external magnetic field (H) in the finite-?/ 
impurity Anderson model, employing the slave boson technique [|TjJ which is generalized to 



the finite-?/ case ||11|| . To our knowledge, there has been no previous study of considering both 
the finite-?/ and the magnetic field effect, simultaneously. We obtain the Kondo temperature 
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and the impurity magnetization as a function of U and H. The impurity magnetization at 
a finite temperature is also determined. In the finite-?/ set of coupled integral 

equations appears even in the lowest order 1/Nf expansion. In order to solve the coupled 
integral equations, we have exploited simplifying approximations which reduce substantial 
computational efforts. This treatment allows one to study easily even anisotropic systems, 
such as those with an external magnetic field or with crystalline fields. 

The paper is organized as follows. In section II, the large-A/ slave boson theory is 
formulated for the finite-?/ impurity Anderson model. In section III, approximation schemes 
used in evaluating the integral equations are described, and the [/-dependence of the Kondo 
temperature is provided. Section IV is devoted to studying the model in the presence of the 
external magnetic field both at T = and at a finite temperature. A special attention is 
paid to the field dependence of the magnetization. Finally, conclusions are summarized in 
section V. 

II. LARGE- A^j SLAVE BOSON TREATMENT 

We consider the finite-?/ impurity Anderson Hamiltonian in a partial wave representation: 

H = Y< £ k c{ m C km + £/ m fmfm + V ( k ) ( C k m fm + f m C km ) + U Y n ™>nm, (1) 
k,m fn k,m m'~>m 

where the indices m,m'(= 1,2, • • ■ Nf) denote spin and orbital quantum numbers. The 
operator c\ m (ck m ) creates (annihilates) a conduction electron in the state (k,m) with an 
energy dispersion of = k = |k|. The /-electrons of the magnetic impurity have in- 
dependent energy £/ m . The corresponding creation (annihilation) operator and the number 
operator are fmifm) and n m = f^fm, respectively. The matrix element V{k) represents 
the hoping of an electron between the conduction band and the impurity level. We assume 
that the hoping takes place only between the states with the same m, and that the V(k) is 
independent of m. The Hamiltonian of Eq.(|l|) is studied in the subspace of at most double 
occupancy, that is, in the subspace of Z ,/ 1 , and f 2 states. This is a good starting point, 
since U is large enough to rule out configurations with more than double occupancies. 
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Coleman QT6J has introduced a single slave boson b to treat the infinite-?/ model. Then 



standard many-body techniques can be utilized in this approach without losing the strongly 
correlated properties of the model. Coleman's method can be generalized to study the finite- 
U model, by introducing a set of 'heavy bosons' d) mm , = —d m , m WU^MKB. The heavy bosons 



represent N f( N f ^ doubly occupied f 2 states, energies of which are E 2m m' = £/ m + £/ / + U. 
In this scheme, the Hamiltonian becomes 

^ = E £fc C km C km + E £ fm fmfm + E E 2mm' d ] mm , d m m' 
k,m m m>m' 

+ E^)( C L&Vm + h.C.)+ £ (cL/i4m' + h.C.) . (2) 

fc,m k,m^m' 

This Hamiltonian commutes with the "charge" operator Q defined by 

Q = tfb + J2 fLfm + E dlm'dmm'. (3) 

m>m' 

Physical quantities of the impurity Anderson model should be obtained with the constraint 
Q — 1. One way to impose the constraint is to add a 'chemical potential' term XQ to ?i, and 
to project it onto the physical subspace Q = 1 by taking A — > oo at the end of calculation 

With this prescription, the partition function and other quantities are expressed in terms 
of the spectral functions p a (uj), a = b,f m ,d mm i corresponding to the J ,/ 1 ,/ 2 states, re- 
spectively. The total partition function can be decoupled to a product Z c Zf, where Z c is 
the partition function of the non-interacting conduction electrons and Zf is the impurity 
contribution given by 

Zf = f°° dw e-*> \p b {u>) + E Pf m (w) + E AW M ] • ( 4 ) 

J_00 \ m m>m' / 

Each of the spectral function is related to the imaginary part of the corresponding retarded 
Green's function G a {uj) : 

1 

p a (oj) = ImG a {u). (5) 

7T 

The Green's functions are written in terms of their self-energies IT, S f m , D mm i : 



G b (co) = (u -n(w) + iS)' 1 (6a) 
GfM = (w - e/m - £ /r » + z^)- 1 (6b) 
Gj mm , (w) = (a> - £ 2 mm' - D mm i(u) + it!)) -1 , (6c) 

with a positive infinitesimal number 5. The self-energies can be evaluated by using the "Non- 
Crossing Approximation (NCA)" which is generalized for the finite- £7 [|lOHl4j1. Recall that 
the NCA has been very successful for the infinite-?/ model ||5| Jl6|Jl7| . The "generalized NCA" 
scheme leads to coupled integral equations containing complicated vertex corrections. For 
this reason, it is practically impossible to solve the equations exactly, and so some simplifying 
approximations are necessary to resolve the problem. 

Instead of using the generalized NCA, we adopt here a simple picture of considering 
only the lowest order diagrams in 1/Nj, that is (1/Nf)°. In fact, crucial finite-?/ effects 
are incorporated already in the lowest order diagrams, as shown below. Figure [j] shows the 
lowest order diagrams for the self-energies of the empty boson and the pseudo-fermion /. 
The self-energy of the heavy <i-boson has no contribution in the lowest order. Diagrams 
in Fig.|l] are calculated by using the standard Feynman rules with the projection procedure 
mentioned above: 

II(u;) = - E / B def{e)G fm {u + e)T m (u;e) (7a) 
vr - J-B 



r m (u J ;e) = l + - £ I* de'f(6')G fm ,^ + 6')G dmml (cu + e + e')r m ,(iu;e'), (7b) 

7T tit \ J — B 

rn' ( rtm I 



with the vertex function, 

E ' 

and 

s/ m (w) = ^ E [ B J £ 'f(e) G <i mm >("+ £ )> ( 7c ) 

D mm ,(u) = 0. (7d) 

Here f(e) is the Fermi-Dirac distribution function, and A = 7rp(e)V 2 (e) is the hopping rate 

of an electron between conduction band and impurity state in a magnetic channel. Here 



p(e) is the density of states of the conduction band. The hopping rate is assumed to be a 
constant A for —B<e<B and otherwise. Eqs.([Ta|) - (|7d|) with Eq.(|6|) complete the 
theory of the lowest order in 1/Nf. The present treatment reproduces at T = the exact 
Nf — > oo results of the variational approach P/Hl- 

In the spin fluctuation limit, where f° and f 2 configurations are energetically unfavorable 
relative to f 1 , most important properties of the Anderson model are determined by the b- 
boson self-energy n. The II contains diagrams representing successive "spin-flip" scatterings, 
which consist of two kinds of processes as shown in Fig|J. These are the processes which 



\v\ IVI 

give rise to the coupling constant J ~ — (ej+u) De t ween conduction electron spin and the 



impurity spin, in the Schrieffer- Wolff transformation | 12| . While the process (a) in the figure 
is present in the infinite-C/ model having an f° intermediate state, the process (b) which 
contains the <i-boson line as an intermediate state does not exist in the infinite-C/ model. The 
process (b), which comes out from the vertex function T m in the empty boson self-energy, 
leads to a renormalization of the Kondo temperature through the modification of the 
coupling constant J(U) and Ta(U) ~ B exp [l/NfJ(U)p(0)]. 



III. APPROXIMATIONS AND THE KONDO TEMPERATURE 

Thermodynamic properties of the finite-?/ impurity Anderson model in the restricted 
subspaces f°, f 1 , and f 2 are determined entirely by the spectral functions p a (u). The ground 
state energy E is obtained from the lowest pole of p&(a>) at T = 0, which corresponds to a 
non- magnetic ground state. Using Eq.([Ta|), one gets the following equation for E : 

= n(E ), (8) 

with II calculated at T = 0, where E is real. Similarly, the lowest poles of Pf m (uj) yield 
magnetic excited energies 8f m satisfying 

e /m -e /m -E /m (e /m ) = 0. (9) 
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S/ m (e/ m ) has no imaginary part in the temperature range of interest, T < Ta- The "Kondo 
temperature" T4, which is defined by the difference between the ground state energy E 
and the lowest excited state energy min{£/ m }, characterizes the low temperature and low 
energy properties of the model. Eq and £/ m calculated in this fashion reproduce the exact 
Nf —> 00 ground state and the excited state energy, respectively, obtained by the variational 
approach [pj|l4|. 

To get Eq, one should solve numerically the set of coupled integral equations involving 
the vertex function T m . Note that in the infinite-?/ case, T m = 1. To solve the equation, we 
have devised two simple approximations, which can be readily applied to anisotropic systems 
without much computational efforts. In T m , there are two independent energy variables, u 
and e. While T m has a large dependence on uj near uj = E , due to Gf ,(uj + e') in the 
integrand of the Eq.([7b]), it has relatively a weak dependence on e. The e dependence of Y m 
exists only in Gd , (uj + £ + a') = (uj + e + e' — E2mm')~ l ■ This dependence would be weaker 
for larger U. Thus we replace T m i(oj;e') in the integrand by r m /(a>;e) to avoid solving the 
coupled integral equations. Actually, this corresponds to neglecting some energy conserva- 
tion of diagrams for vertex function in Fig.l. However, owing to the weak dependence on e 
of the Eq. (|7E|) , this approximation turns out to be quite good, as will be shown below. The 
resulting equation reads 

T m (uj;e) = 1 + — V A mm ,(uj;e)T m ,(uj;e) } (10a) 

where 

A mm ,(u;e) = [ B de' f(e')G fm ,(uJ + e')G dmm ,(uJ + e + e'). (10b) 

Now Eq. (|10a|) is merely a linear algebraic equation of T m 's for given energy variables uj,e. 
This is our first approximation (I). 

The equation can be further simplified by representing the pseudo-fermion Green's func- 
tion as 



where Zf m is the renormalization coefficient at u = if m defined by 



d 



(lib) 

This representation is quite reasonable, since the incoherent background of Gf m for large 
U is located at very high energy (lo > E 2mm ') as compared to if m , and the weight of the 
incoherent part (1 — z/ m ) would be very small (note that it becomes zero for infinite U). 



Thus representing Gf m by Eq.(lla) makes no practical difference in evaluating Eq. ( |10b|) for 
large U. Then A mm / in Eq.(|10b|) can be obtained analytically at T = 0, 



A-kt)^ log ( !%Z^_gw ~"- £ R ). (12) 

This is our second approximation (II). 

The present approximation schemes yield fairly good results of the ground state energy 
Eq and the Kondo temperature T4 for [/ > 5 eV. As shown in Fig|| in the large U limit 
(U B), the Kondo temperature as a function oil/U increases or decreases depending on 
the parameters, £f, A = A^jA, and S 0. Results obtained by using both approximations 
(I and II) are close to those obtained by using the approximation (I) only. Furthermore, 
both results are in good agreement with those of the exact Nf — ► 00 variational approach. 
This demonstration indicates that the approximations (I) and (II) we employed are very 
reasonable at least relatively large U. Deviations from the variational calculations for smaller 
U appear because the approximation schemes made above become less exact for smaller 
value of U, as pointed out above. We used in Fig.|| quite a large value of the half-bandwidth 
B = 24 eV. One can conjecture that the approximation (I) works better for a realistic 
smaller bandwidth, e.g. 2B < O(10eV), because, for smaller B, the e'-dependence of the 
vertex function T m /(uj,e') in the integrand of Eq. ([7b]) would be weaker. 



IV. SCALINGS IN THE PRESENCE OF THE EXTERNAL MAGNETIC FIELD 

With the approximation schemes described in the previous section, we investigate effects 
of the external magnetic field in the finite-?/ impurity Anderson model. By applying the 
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external field H = Hz to a degenerate system (ej m = Ef), Zeeman splitting occurs in the 
localized /-levels; £/ m — £f + gfi B Hm, m = —j, — j + 1, • • • , j, (2j + 1 = Nf). Conduction 



electron polarization is neglected, since its contribution is perturbatively small |18[ . 

The Kondo temperature in the presence of H is also obtained from the energy difference 
between the ground state and the lowest excited state energy. Figure (| plots the magnetic 
field dependence of Tjy. TVs are normalized by their zero-field values as given in Table I. It 
is seen that T4 decreases with increasing H. A reduction of the Kondo temperature by the 
applied field originates from the reduction of the ground state binding energy. The applied 
field lifts up the degeneracy of the /-level to decrease the effective degeneracy. Hence the 
binding energy will be reduced exponentially, recalling that T4 ~ B exp in the Kondo 
limit of the infinite-f/ model. Notable in this figure is that T a (H)/Ta(0) exhibits almost the 
same functional form against gusH/ksTAiO) regardless of the size of U, implying that the 
scaling behavior holds with high accuracy in the low field region. In contrast, Ta{H)/Ta{0) 
has a different scaling behavior for different Nf. In the infinite-f/ limit, Ta(H)/Ta(0) can 
be expressed as (see Appendix), 

T A (H) 1 



T A (0) fff-i 

1 lm=l 



1 , m g^ B H/T A (0) 
L "T" " L T A (H)/T A (0) 



l/N f 



exp 



l(E o (H)-E o (0)) 



(13) 



The solution of Eq.([T3|) can be regarded as an ^/-dependent scaling function characterizing 
the low T and low H properties of the model. Using the Kondo temperature T A rescaled by 
U, the scaling of Eq. (|T3"| ) will be also valid in the finite-?/ case. 

The impurity contribution to the magnetization M imp of the system can be obtained 
from the relation 

M ' mP = ^ l0gZ/ - (14) 
The ground state is no longer non-magnetic in the presence of the magnetic field, since 
the field polarizes both the impurity electrons and the conduction electrons. In Fig.|5|, the 
ground state impurity magnetization M zmp (H) is plotted as a function of q^bH ^//cbT^(O). 
The magnetization in this figure is normalized by its saturation value M*^ p = jg^iB- The 



magnetization at low field also shows a scaling behavior almost perfectly, as in the case of 
Ta{H). The present results of M imp (H) are compared with the Bethe ansatz results for 
the Coqblin-Schrieffer model JO?) . In the Bethe ansatz results, appears the low temperature 
scale Tl which is related to the zero temperature susceptibility xo = {g^B) 2 j{j + 1)/3Tl- 
On the other hand, the susceptibility in our treatment is given by Xo = n f ^ 9(J,B ^J^ +1 ^ for 
U = oo. By considering the relationship between two energy scales, it is possible to compare 
our results with the Bethe ansatz results. As is seen in Fig.[5], the present results of M imp (H) 
agree quite well with the Bethe ansatz results in the low field region, even though they are a 
bit overestimated in the high field region. The deviation from the exact results at high field 
indicates that the 1/Nf expansion may not work well in the high field region. As mentioned 
previously, the applied field splits the /-level and reduces the effective degeneracy, and so 
the lowest order 1/Nf expansion becomes invalid in the high field region. It is noteworthy 
that the overestimation of M imp (H) observed in Fig.|5| is also found in the infinite-f/ NCA 
scheme for large H/T pOf . 

An interesting feature is found in the low field magnetization. The inset of Fig.^ indicates 
that the low field magnetization increases more rapidly than linearly in H. This deviation 
from the linearity, so called the "superlinear" behavior, has been reported by Hewson et 



al ||19|| , based on the Bethe ansatz solution of the Coqblin-Schrieffer model. They found 
that the superlinearity occurs for Nf > 3, and that the experimental results of YbCuAl 
are well described by their Nf = 8 model. The superlinearity is also revealed by the finite 



temperature NCA scheme as well as by the mean field approximation |2T[ in the infinite- 
U, large- iVy treatments of the Anderson model. The superlinear behavior has its origin in 
the fact that the Kondo resonance, which is located near Ta above the Fermi level, becomes 
narrower as Nf increases. That is, the Kondo resonance is sharply defined for large Nf, 
implying a strong screening effect. As a consequence, the spin polarization does not begin 
to dominate, until the magnetic field becomes comparable to Ta- As the Zeeman energy is 
comparable to I^bTa, the magnetic field overcomes the screening by conduction electrons, 
and the spin polarization dominates the system. It is quite natural that the superlinearity 
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is present in our finite-?/ model, in view of the almost perfect scaling of M imp (H) in the low 
field regime. The super linearity persists at low T/T4, but disappears at higher temperature 
because of suppression of the Kondo resonance. 

At finite temperature, the situation is more intriguing. For a given low T, M imp (H) shows 
a reupturn at a certain field H* (denoted by arrows in Fig|]), apart from the superlinearity 
at lower field. Making more investigation into this anomaly, one could find that the Kondo 
temperature at the reupturn point is nearly equal to a given temperature of the system, 
that is, Ta{H*) ~ T. This phenomenon could be understood as follows. For a given low T, 
M mp (H) has a superlinearity at low H, but diminishes its increasing rate for intermediate 
values of H. As H increases further, M tmp (H) exhibits a reupturn at H = H* which fulfills 
Ta(H*) ~ T. It takes place because the population in the lowest excited magnetic state 
increases abruptly at the point satisfying Ta(H*) ~ T, producing observed reupturns at H*. 
It is shown in Fig.|6| that the reupturn smears out at higher temperature, since the Kondo 
resonance is suppressed at high temperature. Another thing to note is that the scaling of 
M imp (H) at finite T is not so good as that at T = 0. 

Figure |7| present M tmp (H) for U = 5eV at various temperatures. Non-monotonicity as a 
function of T is apparent in the low field magnetization. With increasing the temperature, 
M mp (H) for a given small H, which is proportional to magnetic susceptibility, increases 
first at low T, but after reaching a maximum at an intermediate T*, M imp (H) decreases 
at higher T (see the inset of the Fig.0). This non-monotonic behavior in the temperature 
dependence of the magnetic susceptibility is known to be more pronounced for a system 
with larger Nf |22], and is considered to have the same origin as that of the superlinearity 



in M imp (H) 0. 



V. CONCLUSION 



We have studied the finite-C/ impurity Anderson model in the presence of an external 
magnetic field. Large- Nf expansion considering f°, f 1 , and f 2 subspaces in the slave boson 
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representation is employed to take into account finite-?/ effects. The lowest order treatment 
gives rise to a complicated vertex function in the empty boson self-energy and a renormal- 
ization of /-level energies. We have devised two simple approximation schemes to evaluate 
the vertex function. One consists in transforming the integral equation into an algebraic 
equation by assuming that the vertex function has a weak dependence on one of two energy 
variables. The other is to simplify the pseudo-fermion Green's function by neglecting the 
incoherent part. These approximations reduce computational efforts substantially, and are 
shown to yield fairly good results for U ^ 5 eV. 

This treatment has been applied to the system with an external magnetic field. It is 
found that the Kondo temperature is reduced by the magnetic field. Scaling behaviors in 
the field dependence of the Kondo temperature and the impurity magnetization are found to 
hold almost perfectly at low field and low temperature. This implies that the main effect of 
the finite-?/ in the presence of the magnetic field is manifested through the renormalization 
of the Kondo temperature. Some intriguing features are found in the field dependence of the 
magnetization, such as superlinearity at low field and reupturns at a higher field H* for a 
given low temperature, which are expected to occur from a competition between the singlet 
binding energy and the Zeeman energy gain of electrons. 
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APPENDIX: 

Let's consider the Kondo temperature in the presence of the applied magnetic field for 
U — oo case. Degenerate /-levels split by the magnetic field. The energy levels can be 
expressed in terms of the Kondo temperature, as defined in section III : 

12 



£ fm = Eo(H) + T A (H) + m gf i B H, (Al) 
where m = 0, 1, 2, ■ • • , Nf — 1. In the U — > oo limit, Eq.(^) reduces to 

ft< g >~7g/ U( g) - E - £/ . < A2 > 

at T = 0. Analytical evaluation of Eg. ( |A2[ ) with Eg. (|AT|) yields 

«m-f^{'ff(i + ^)p. 

In the derivation, we used the condition, B ^> Nfgfi B H and B ^> T4. Eg . ( |A3|) corresponds 
to Eg.©. 
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FIGURES 

FIG. 1. Lowest order ((1/iV/) ) diagrams of the empty boson and the pseudo-fermion 
self-energies in the finite-?/ Anderson model. The full, wiggly, dashed, and jagged lines stand 
for conduction electron, empty state boson, pseudo-fermion, and heavy boson propagators, respec- 
tively. Double dashed line denotes the renormalized pseudo-fermion propagator as represented by 
the diagram in the bottom. 

FIG. 2. Two elementary "spin flip" scattering processes appearing in the empty boson 
self-energy II, having an ^-intermediate state (a) and an / 2 -intermediate state (b). 

FIG. 3. Comparison of the present results for the Kondo temperature with the exact Nf — > oo 
variational results at T = (Ref. [9]). Dashed (solid) lines are the Kondo temperatures obtained 
by the approximation scheme (I)((I) and (II)), and points correspond to the exact variational ones. 
Parameters used here are Nf = oo, A = NfA = 0.75 eV, B = 24 eV with £/ = — 2.5eV in (a) and 
£f = — 1.5eV in (b). Here degenerate /-level Ef m = £f is assumed. 

FIG. 4. Normalized Kondo temperature vs. normalized magnetic field strength for Nf = 4 
and Nf = 6. The normalized Kondo temperatures are almost the same for different f/'s. Other 
parameters used are £/ = — 2eV, A = 0.9 eV, B = 3eV. 

FIG. 5. Scaling behavior of the impurity magnetization which is normalized by its saturation 
value M^ t p = jg/J-B- Bethe ansatz results of the Coqblin-Schrieffer model for Nf = 4(+) and 
Nf = 6(0) are also given, for the comparison. The inset shows deviations of the low-field magne- 
tizations from linearity for Nf = 4 and Nf = 6. Solid lines are the present results for U = 5eV, 
and dashed lines are extrapolated values from the initial slope, ^ h. Other parameters are 
the same as those in Fig.4. 
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FIG. 6. Field dependence of the magnetization for N f = 4 at (a) T/T A = 0.1, (b) T/T A = 0.3, 
(c) T/Ta = 0.5, and (d) T/Ta = 0.9. Other parameters are the same as those in Fig.4. Arrows 
in (a), (b), and (c) denote reupturns of the magnetization. The reupturn smears out at higher 
temperature as shown in (d). 

FIG. 7. Field dependence of the magnetization for various temperatures with Nf = 4 
and U = 5eV. Other parameters used are the same as those in Fig.4. The inset shows the 
non-monotonic behavior of the magnetic susceptibility divided by its value of T = 0, x(T)/x(0). 
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TABLES 

TABLE I. Kondo temperature Ta at H = for Nf = 4 and Nf = 6. Parameter sets are the 
same as those in Fig. 4. 

U (eV) T A (meV) (N f = 4) T A (meV) (N f = 6) 

oo 2.76 2.76 

10 6.04 6.51 

5 11.9 13.4 
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